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Система мониторинга и морфология 
формы нестационарных сигналов 
внутричерепного давления 


В статье рассматривается разработанная и внедренная авторами система измерения внутричерепного 
давления (ВЧД). Предлагается для изучения волновых форм временных сигналов ВЧД использовать 
алгоритмы калмановской фильтрации и калмановского сглаживания, которые в каждый момент времени 
адаптивно производят оценку параметров авторегрессионной модели. С помощью предложенных 
алгоритмов изучена морфология волн ВЧД при острой травме головного мозга. 


Введение 


Измерение внутричерепного давления (ВЧД) является крайне важным при ней- 
рохирургическом лечении. По оценкам Американской ассоциации нейрохирургов [1] 
только в течение 2002 года общее количество измерений ВЧД составило 110 000 
процедур. Внутричерепная гипертензия наблюдалась в более чем 50% пациентов с ост- 
рой травмой головного мозга (ТГМ) [2]. Кроме того, частота случаев повышенного ВЧД 
в течение госпитализации у пациентов располагалась в диапазоне от 30% до 80% [2]. 
Использование мониторинга ВЧД позволяет получить наиболее вероятный прогноз раз- 
вития патологии, является методом, который позволяет провести патофизиологически 
обоснованный комплекс лечебных мероприятий [3]. 

Данные мониторинга ВЧД определяют методы снижения давления, например, пу- 
тем дренирования спинномозговой жидкости (СМЖ), что приводит к улучшению моз- 
говой перфузии. 

В мировой практике разработано и внедрено множество способов контроля ВЧД. 
Системы мониторинга ВЧД классифицируются по месторасположению или по типу сис- 
темы датчика преобразования давления. Фондом травмы головного мозга (Вгаш Тгаита 
Еоипданоп) совместно с Американской ассоциацией нейрохирургов (Атепсап Аззос1айоп 
ог Меиго|о?са| Зигоеоп$) был опубликован обзор работ по технологии мониторинга 
ВЧД [4]. Оптимальным внутричерепным местоположением систем мониторинга ВЧД 
являются вентрикулярное, паренхимальное, субдуральное, субарахноидальное и эпиду- 
ральное пространства. Вентрикулярное ВЧД используется как исходный эталон для 
сравнения точности измерения давления в других внутричерепных отделах. 

Ранжирование систем мониторинга ВЧД, основываясь на их точности, надежно- 
сти и стоимости, имеет следующий порядок [4]: интравентрикулярные системы — жид- 
костно-заполненный катетер с внешним тензометром; интравентрикулярные системы — 
микрокристаллический тензометр или волоконно-оптический датчик; паренхимальные 
системы преобразования давления; субдуральные системы; субарахноидальные жид- 
костно-заполненные системы и эпидуральные системы. После имплантации катетера 
или датчика система передает изменения давления к измерительному преобразователю, 
который конвертирует импульсы в видимую индикацию, или в цифровые величины. 
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Целью данной работы является разработка и исследование предлагаемой авто- 
рами наиболее точной по измерению и дешевой по стоимости, вентрикулярной жидкост- 
но-заполненной системы для мониторинга ВЧД с внешним тензометром, позволяющей 
производить калибровку ш у1уо. Использование алгоритмов калмановской фильтрации 
и калмановского сглаживания, которые в каждый момент времени адаптивно произво- 
дят оценку параметров авторегрессионной модели. С помощью предложенных алгорит- 
мов изучена морфология волн ВЧД при острой травме головного мозга (ОТГМ). 


Материал и методы исследования 


Проведено обследование 19 пострадавших с тяжелой черепно-мозговой травмой 
(глубина коматозного состояния по шкале ком Глазго менее 8 баллов), находящихся на 
лечении в отделении нейрореанимации Донецкого областного клинического террито- 
риального медицинского объединения. 


Внешний жидкостно-заполненный тензометрический 
(пневматический) измерительный преобразователь 


Наиболее распространенная жидкостная система мониторинга ВЧД представляет 
собой катетер с эластичным баллончиком на ее конце, полностью заполненным жидкой 
средой, которая передает давление от области измерения на мембрану преобразователя 
давления. При импульсном воздействии на жидкостную замкнутую систему катетер-пре- 
образователь могут возникать собственные колебания. Частота и коэффициент демпфиро- 
вания таких колебаний определяется тремя механическими коэффициентами: податли- 
востью (С), инерцией жидкости (1) и сопротивлением жидкости (К). Податливость свя- 
зана с гибкостью мембраны преобразователя давления и эластичностью баллончика 
катетера. Инерция жидкости определяется коэффициентом, обусловленным давлением 
ускорения жидкости через катетер. Сопротивление жидкости зависит от вязкости жид- 
кости в катетере и ее трения о стенки катетера. Резонансная частота Р, и коэффициент 
демпфирования & имеют следующий вид: 


о 
ке: 4-9. С 


Из этих формул следует, что при малых коэффициентах податливости и инерции 
жидкости резонансная частота сдвигается в сторону высоких частот. 

Полученный коэффициент демпфирования условно можно разбить на три части: 
слабо демпфированный, сильно демпфированный и оптимально демпфированный. В сла- 
бо демпфированной системе волны давления имеют тенденцию отражаться в пределах 
катетера и трубки преобразователя, создавая гармоники стоячих волн. В сильно демп- 
фированной системе фрикционные силы препятствуют передаче волны ВЧД таким об- 
разом, что она теряет энергию. В слабо демпфированной системе из-за сдвига фаз между 
колебаниями ВЧД и собственными колебаниями (стоячими волнами) форма волны ВЧД 
становится узкой, достигая пикового максимума, что приводит к ошибкам до 50% 
в измерении ВЧД. С другой стороны, с ростом демпфирования уширяется и сглаживается 
форма волны ВЧД. 

Физические свойства катетера также влияют на передачу формы волны давления. 
Как правило, более жесткий, короткий и белее тонкий катетер приводит к передаче более 
точной формы волны давления [5]. 

Мы изготовили жидкостную систему катетер-преобразователь давления (тензометр), 
где полиуретановый катетер с более эластичным силиконовым баллончиком на конце за- 
полнен физиологическим раствором и присоединен к тензометру через капиллярный блок. 
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Капиллярный блок позволяет сепарировать колебания внутричерепного давления 
от собственных колебаний жидкой среды в тензометре. Заполненный физиологическим 
раствором полиуретановый катетер длиной 250 мм и внутренним диаметром 1 мм с си- 
ликоновым баллончиком на конце является оптимально демпфируемой системой. В тен- 
зометре средой, передающей давление, является вязкая силиконовая жидкость. Система 
капиллярный блок-мембрана тензометра представляет собой избыточно демпфируемою 
систему, в которой подавляются собственные колебания, частота которых, согласно (1), 
сдвигается в сторону высоких частот и коэффициент затухания растет из-за увеличения 
вязкости силиконовой жидкости. Более того, капиллярный блок позволяет сепарировать 
стерильный физиологический раствор в катетере от вязкой силиконовой жидкости в 
тензометре, не допускает их смешивание. 

В капиллярном блоке предусмотрен узел, позволяющий производить калибровку 
датчика ш У1уо. Контрольной точкой является атмосферное давление. Обнуление про- 
исходит открытием запорного крана на атмосферное давление. Запорный кран присо- 
единяет к капиллярному блоку либо катетер, либо воздушный интерфейс. Обнуление 
устраняет нулевой дрейф, который возникает при смещении катетера в результате дви- 
жения пациента. 

Измерение давления производится с помощью сенсора, принцип действия кото- 
рого основан на тензорезистивном эффекте — изменении электрического сопротивления 
полупроводника в результате действия нагрузки, создающей деформацию. Чувстви- 
тельным элементом сенсора является мембрана с тензорезисторами, сформированными 
в эпитаксиальной пленке кремния на подложке из сапфира, припаянной серебросодер- 
жащим припоем к кварцевой мембране. Терморезисторы соединены в мостовую схему. 
Под действием давления измеряемой среды мембрана прогибается, тензорезисторы ме- 
няют свое сопротивление, что приводит к разбалансировке моста Уинстона. Разбаланси- 
ровка линейно зависит от степени деформации резисторов и, следовательно, от приложен- 
ного давления. Мы использовали тензорезистивный сенсор разработки Института физики 
полупроводников НАН Украины, НПФ «Рост», г. Киев. Для уменьшения температурной 
погрешности измерения давления была предусмотрена пассивная схема термокомпен- 
сации, состоящая из нескольких напыленных тонкопленочных компенсационных термо- 
резисторов, снижающих температурный дрейф нуля выходного напряжения в 6 - 10 раз. 

На основе жидкостно-заполненной системы катетер-тензометр нами был разра- 
ботан прибор ВЧД мониторинга с измерительным модулем М3500Р (УкКрРЦЕНДИСИ 
НАН Украины, г. Киев) [6]. В настоящее время прибор используется в нейрохирургии 
для мониторинга нарастающей внутричерепной гипертензии у пациентов с черепно-моз- 
говой травмой, ишемическим и геморрагическим инсультом, опухолью мозга [7], [8]. Си- 
стемное программное обеспечение позволяет адресное обращение к измерительным моду- 
лям, введение информации к компьютеру, накопление и математическую обработку дан- 
ных и их графическое изображение в течение мониторинга несколько суток непрерывно. 


Адаптивная авторегрессионная (ААР) модель 
нестационарных колебаний ВЧД 


Наблюдаемые сигналы внутричерепного давления представляют собой временной 
ряд, который математически может быть представлен авторегрессионной (АР) моделью. 
Естественные сигналы (речь, музыка или сейсмические сигналы), а также многие био- 
медицинские сигналы (электроэнцефалограмма) имеют авторегрессионную структуру. 
Однако эти сигналы искажены совокупным широкополосным шумом (шумом измерения, 
экологическим шумом), следовательно, необходимо усиление сигналов для дальнейшей 
их обработки. 
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Авторегрессионная модель выводит путь, по которому синтезируется вероятност- 
ный процесс временного ряда х(и). Этот процесс х(и) можно создать, применяя меро- 
морфный (все полюса) фильтр к белому шуму е(п) с нулевым средним и дисперсией о”. 
Этот процесс схематично показан на рис. 1. 


Рисунок 1 — АР процесс, где е(п) является входом (нулевое среднее белого шума) 
и х(п) является выходом (сигнал, который будет смоделирован) 


Следует обратить внимание, что сигнал при анализе является фактически выходом 
системы. Это является существенным отличием АР модели от других методов спект- 
рального анализа, в которых исследуемый сигнал разлагают на составляющие при входе к 
фильтру (например, быстрое преобразование Фурье). Поэтому коэффициенты фильтра 
содержат информацию относительно спектрального содержания х(и) и определяются как 
параметры АР модели. В соответствии с этим АР модель является параметрической мо- 
делью. 

Модель АР весьма проста и полезна для описания стохастического поведения 
временного ряда. Он описан следующим уравнением [9]: 


Р 
х(п) = У`арх(и-К)+е(п), (2) 


К=1 
где х(и)— анализируемый физиологический сигнал (сигнал ВЧД) в момент и, {а и — 


параметры модели, ИС _ ЮР — задержанные выборки сигнала и е(п) — вход модели, 


который принимаем как независимую случайную последовательность и обычно 
распределенную с нулевым средним (белый шум) и дисперсией о”. Индекс и является 
целым числом и описывает дискретные, равноудаленные моменты времени. Время # в 
секундах является {Е = и//, =ихАТ с частотой выборки / и интервалом выборки ДТ = 1/№. 
Для простоты вводится векторная запись для АР параметров и вектора, Х(и)”: 


а(п) = [@,(п),..а, т, 
Х(иТ = [(®),....жи-р+ ОГ, (3) 


где а/(И), ..., ар(п) — авторегрессионные коэффициенты при времени и; х(и-р-+/), ..., х(п) — 
последовательные (поочередные) данные от анализируемого сигнала ВЧД. Соответ- 
ственно модель АВ может быть написана как 


х(п) = а(п)" Х(п)" + е(п) (4) 
и передаточная функция в /-области имеет вид 
х(2)/е(2) =/(-ах2'-...-а,х2?). (5) 


Спектральная плотность х(п), обозначенная как 5лк, связана с АР параметрами 
для модели р-го порядка уравнением [10]: 


от 


2 
Р —72МАТ 
|- > + ак (пе рт 


блв(и, Л) = (6) 
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Авторегрессионный метод анализа является параметрическим методом, обладаю- 
щим большой вычислительной мощностью. Поэтому в состоянии описать полностью 
статистику второго порядка временного ряда, включая спектральный анализ. АР-спектр 
является спектром максимальной энтропии и, следовательно, описывает функцию 
спектральной плотности наиболее точно с минимальным числом параметров. Кроме 
того, АР моделирование не требует выбора определенного диапазона частот. Индиви- 
дуальные изменения средней частоты, которые могут быть важными при исследовании 
сигнала, можно более легко рассмотреть в АР модели. 

Обычно ВЧД является неустановившимся процессом, и АР параметры изменяются 
со временем (так же, как и дисперсия оз). Для адресации нестационарного поведения 
биомедицинских сигналов необходимо использовать окно смещения оценок АР пара- 
метров. В зависимости от размера и сдвига окна традиционные алгоритмы оценки АР 
параметров могут требовать болыших вычислительных затрат. Вместо этого можно оце- 
нить нестационарные АР параметры, вводя динамическую модель и используя фильтр 
Калмана [11]. Такой подход называется как адаптивное авторегрессионное (ААР) моде- 
лирование и успешно использовался при анализе ЭЭГ [9]. 


Калмановская фильтрация нестационарных колебаний ВЧД 


Фильтр Калмана был предложен в 1960 году доктором В.Е. Ка|тап, как оптималь- 
ный рекурсивный фильтр для оценки вектора состояния из данных измерения искажён- 
ных шумом [12]. Его преимущества перед другими фильтрами, например, такими, как 
фильтр Винера в уменьшении математической сложности обработки больших последо- 
вательностей данных. Блок-схема фильтра Калмана показана на рис. 2, где ОШЗ пред- 
ставляет одношаговую дискретную задержку. Дискретно-временной фильтр Калмана 
обеспечивает рекурсивное решение замкнутой формы для оценки линейных дискретно- 
временных динамических систем, которые могут быть описаны уравнениями вида [10]: 


аи = А Я@т-1 + Чи-—1 
а Вы: (7) 


где а, — обновленный вектор состояния системы на временном шаге и; а». — вектор 
состояния системы на временном шаге и-/; х„ — вектор измерения на временном 
шаге п; О». > №0, О.п — шум процесса на временном шаге п-; у„ = №0, Г) - шум 
измерения на временном шаге и; А»„.; — матрица перехода динамической системы; А’, — 
матрица модели измерения. 

Предшествующее (априорное) распределение для состояния - а - Мао, Ро), где 
параметры ао и Ро устанавливаются, используя информацию, известную о системе 
при исследовании. 


&(пп-1) а(п-1п-1) 


Рисунок 2 — Блок-схема фильтра Калмана 
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Фильтр Калмана функционирует в два шага: шага предсказания, где следующее 
состояние системы предсказано беря предыдущие измерения, и шага обновления, 
где текущее состояние системы оценено беря измерение при этом временном шаге. 
Шаги приводят к следующим уравнениям: 

Шаг предсказания: 


ал = и Я, (8) 
те Т 
В = и—1 Ру О 1. (9) 
Шаг обновления: 
ве (10) 
е„ — разность между фактическим и расчётным измерениями на временном шаге и. 
В, =Х,Р; Х/ +И,, (11) 
К„ - ковариация предсказания измерения на временном шаге и. 
—1 
К, =Р.Х, В” (12) 
К,- коэффициент усиления фильтра. 
ар =а, — Куре (13) 
ВЕР -К В КГ (14) 


предполагаемое среднее а„ и ковариация состояния Р, на временном шаге и после 
наблюдения измерения. 

Эти уравнения используются рекурсивно для получения оценки состояния при 
данном и. Существует два основных фактора, которые могут влиять на работу фильтра 
Калмана [10], первый — введение параметров фильтра Калмана, таких, как ковариация 
шума процесса, ковариация шума измерения, и второй — начальные условия. Эти па- 
раметры — настраиваемые механизмы фильтра. 

Из приведенных уравнений фильтра Калмана видно, что усиление зависит от 
ковариации предсказания и ковариации шума измерения И». Ковариация предсказания 
также зависит от ковариации шума процесса О„. Если ковариация шума процесса в урав- 
нении (9) большая, увеличивается ковариация предсказания. Поэтому усиление фильтра 
Калмана в уравнении (12) можно рассматривать как баланс между ковариацией шума 
процесса к ковариации шума измерения. Учитывая это обстоятельство, с ростом кова- 
риации шума процесса возрастает усиление фильтра Калмана и, следовательно, увели- 
чивается полоса пропускания. Это вызывает более быстрый переходный отклик, кото- 
рый приводит к большему количеству шума в оценках, произведенных уравнением (8). 
Такой же эффект может быть достигнут, уменьшая ковариацию шума измерения в урав- 
нении (11). Если ковариация шума процесса будет минимальной, произойдет противо- 
положный эффект, который подразумевает, что меньше шума будет присутствовать 
в оцененном состоянии. 

Что касается выбранных начальных условий, только часть действующего алгорит- 
ма будет переходной частью. Поскольку обрабатывается большое количество данных, 
начальные условия исчезают, в конечном счете достигается устойчивое значение сос- 
тояния. Выбирая большую ковариацию предсказания, больше акцента будет уделено 
измерению и меньше - к модели в переходной фазе. 


Оптимальный сглаживающий фильтр с фиксированным интервалом 


Задача сглаживания отличается от фильтрации тем, что для извлечения информа- 
ции о сигнале в момент времени и используются данные, полученные позже. Это озна- 
чает, что в процессе сглаживания в получении результата присутствует задержка. Процесс 
сглаживания в статистическом смысле является более точным, чем фильтрация, посколь- 
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ку при сглаживании можно использовать данные, полученные не только после момента 
времени и, но и до него. Алгоритм оптимального сглаживающего фильтра с фиксирован- 
ным интервалом, известный как сглаживающий фильтр Каисб-Типо-56лебе| или дис- 
кретно-временной сглаживающий фильтр Калмана, был разработан для использования 
в сопряжении с фильтром Калмана [11]. Алгоритм этого фильтра повторно вычисляет 
каждую оценку, произведенную фильтром Калмана, основанную на информации, полу- 
ченной для всего набора проанализированных данных. Эта процедура выполняется блок- 
схемой, показанной на рис. 3, где ДВШ обозначает один дискретно-временной шаг. 
а(1-1п-1) 


а(п-11м> 


а(пп-1) 


- двш 
аа 


Рисунок 3 — Блок-схема оптимального сглаживающего фильтра 
с фиксированным интервалом 


Среднее значение и ковариация а,’ и Р‚° вычисляются с помощью следующих 
уравнений: 


ан = Аха, , (15) 

Р.1 =А„Р, Ал + Ор, (16) 
би =Р, АР, (17) 
ау =а, + [ал +1 —а,+1], (18) 
РР, +5 [Ру — Ра] $, (19) 


где в и Р5 — сглаженные оценки средних значений состояния и ковариации 
состояния на временном шаге и; а, и Р, — оцененные фильтром средние значения 
состояния и ковариации состояния на временном шаге и; а,-1 и Р„+,— предсказанные 
средние значения состояния и ковариации состояния на временном шаге и+/ те же, 
что и в фильтре Калмана; $, — усиление сглаживающего фильтра на временном шаге 
п, которое указывает, насколько сглаженные оценки должны быть подправлены на 
этом определенном временном шаге. 

Из этих уравнений следует, что рекурсия в фильтре Калмана продвигается вперед, 
а в сглаживающем фильтре назад. В сглаживающем фильтре рекурсия начинается с 
последнего временного шага М№ с а" = ами Р8 = Ру. В этом и состоит основное отличие 
между фильтром Калмана и сглаживающим фильтром Калмана. 


Результаты измерений и обсуждение 


Лундберг (Гли4Бег®) идентифицировал три различных типа динамики изменений 
среднего значения ВЧД, как А-, В- и С- волны. Проводился детальный анализ В-волн, 
поскольку их наличие может свидетельствовать об увеличении риска развития внутри- 
черепной гипертензии. Неустойчивый характер этих волн, визуальный анализ временного 
сигнала, искажения волн обычными преобразователями, неадекватное осуществление 
выборки и плохое разрешение монитора часто приводят к исчезновению волны. Исполь- 
зуя вышеописанную методику, мы наблюдали у пациентов с острой травмой головного 
мозга различные типы В-волн — синусоидальной формы (рис. 4), пилообразной формы, 
их суперпозицию с пульсовыми и респираторными волнами (рис. 5 и рис. 6). Наблю- 
даемые нами медленные ВЧД волны попадают в три отчетливых частотных диапазона, 
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определенных в [13], как ультра В-волны от 0,05 Гц до 0,2 Гц; В-волны от 0,008 Гц до 
0,05 Гц; и инфра В-волны ниже 0,008 Гц. 


а) 
Зюпа!+\м/п не Саизэап позе 
21 т 

[= 
== 
Е 
Е 20.5 
9 
= 
2 
= 201 
о. 
= 
= 

49.5 1 : 1 | : 

( 20 40 50 80 100 120 140 160 180 200 220 
Тите (зесопа) 
ЗтооШефд эдпа! 

600 
<> 400 
в-3 
= 
3 200 
| 
В. 

-200, | А Г 

( 20 40 60 80 100, 120 140 160 180 200 220 
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Ка!тап Ез#гтаед Эюпа! Тке ите уагупд зресгит 


0.2 
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> 
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э 
> 
[23] 


“ 100 200 300 80 100 120 140 16С 
Роме\мт Зресга! ВепзКу Тте (5ес) 
Рисунок 4 — Пример медленной синусоидальной В-волны при ТГМ, оцененный 
сглаживающим фильтром Калмана (больной Н-ов): а) ВЧД сигнал; 
Ь) спектральная плотность мощности и изменяющийся во времени спектр В-волны 


В нормальном физиологическом состоянии ВЧД показывает пульсирующую реги- 
страцию на двух различных частотах — одно синхронное с артериальным пульсом, в то 
время как другое медленнее, согласованное во времени с дыханием (респираторной 
компонентой). Обычно амплитуда дыхательных колебаний больше, чем сердечных коле- 
баний, но при внутричерепной гипертензии артериальные пульсации также принимают 
большую амплитуду. Сосудистые волны вызваны артериальными пульсациями в больших 
сосудах внутри головного мозга и производят колебания в объеме желудочковой системы. 
Дыхательная (респираторная) волна синхронна с изменениями в центральном венозном 
давлении и отражает внутригрудное давление. Поэтому В-волны ВЧД почти всегда со- 
держат компоненты и пульсовых волн (рис. 5), и респираторных волн (рис. 6). 

У всех пациентов с ОТГМ наблюдается триггерное переключение волн ВЧД, 
когда волна длительностью 80 — 120 секунд с области частот 0,11 — 0,24 Гц импульсно 
переключается на частоты 0,45 — 0,57 Гц и затем таким же образом переключается в ис- 
ходное состояние (рис. 7). Происходит как бы частотная модуляция волны ВЧД в виде 
меандра с длительностью и амплитудой несколько различной в отдельности для каждого 
пациента. У большинства пациентов мы наблюдали это явление перед возникновением 
ультра В-волн, а иногда и в течение В-волны. 

Достоверная причина такого триггерного переключения волн ВЧД не ясна. Одно 
из возможных предположений состоит в следующем. Наступление волны можно описать 
сосудорасширяющей каскадной моделью, представляющей замкнутый контур нестабиль- 
ного процесса, сопряженного с активным расширением кровеносных сосудов, ростом 
объема мозгового кровотока (ОМК), поднятия ВЧД и уменьшения церебрального пер- 
фузного давления (ЦПД). Относительно быстрое окончание волны может быть описано 
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вазоконстрикцийной каскадной моделью, в которой активное сосудосуживающее собы- 


тие обусловлено уменьшением ОМК и ВЧД и ростом ЦПД [14]. 
а) 
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Рисунок 5 — Пример суперпозиции волн — синусоидальной В-волны и пульсовой 
волны при ТГМ, оцененных сглаживающим фильтром Калмана (больной Н-ов): 
а) ВЧД сигнал; Ъ) спектральная плотность мощности и изменяющийся 
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Рисунок 6 — Пример суперпозиции волн — пилообразной В-волны и респираторной 
волны при ОТГМ, оцененных сглаживающим фильтром Калмана (больной Т.): 
а) ВЧД сигнал; 6) спектральная плотность мощности и изменяющийся 


во времени спектр респираторной и В-волны 
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Рисунок 7 — Пример триггерного переключения респираторных волн ВЧД (больной С.) 


Поскольку сосудорасширяющая/сужающая каскадная модель постулирует интакт- 
ное прежнее регулирование давления до начала волны [14], возможно, что в течение пе- 
риода перед волной критическая частота колебаний передачи цереброваскулярного 
давления будет демпфировать [13]. Однако начало волны, активное расширение кровенос- 
ных сосудов, будет требовать от сосудистого тонуса цереброваскулярного русла резуль- 
тирующего роста в способности цереброваскулярного русла колебаться (вибрировать). 
При окончании волны сосудорасширяющая/сужающая каскадная модель описывает по- 
нижение в ОМК в результате активной вазоконстрикции, как основы для понижения 
ВЧД. Таким образом, тонус цереброваскулярного русла будет возрастать. 


Выводы 


В данной работе представлена разработанная авторами жидкостно-заполненная 
система вентрикулярного мониторинга ВЧД и показано применение фильтра Калмана и 
сглаживающего фильтра Калмана для оценки меняющегося во времени спектра. Оценка 
спектра основана на параметрической ААР модели. В случае, где все наблюдения доступ- 
ны перед применением процедуры оценки, сглаживающий фильтр Калмана может исполь- 
зоваться вместо фильтра Калмана. Однако для обработки в реальном масштабе времени 
использовать эту процедуру нецелесообразно. Показана целесообразность использования 
фильтра Калмана при использовании модели, которая зависит от некоторых параметров. 
Например, в ААР модели допускается изменение параметров во времени, тогда оценка 
этих параметров может быть выполнена, используя фильтр Калмана. Если предположе- 
ние линейности не может быть использовано, в том случае могут использоваться неко- 
торые модификации фильтра Калмана, такие как расширенный фильтр Калмана. 

Необходимы дальнейшие исследования для выяснения роли медленных волн ВЧД 
при мониторинге пациентов с внутричерепной гипертензией, гидроцефалией, тяжелой 
черепно-мозговой травмой, исследования причин обнаруженного триггерного переклю- 
чения волн ВЧД. Перспективно также использовать широкий выбор математических ин- 
струментов для анализа как в частотной, так и временной областях волн ВЧД. 
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Д.В. Дроботько, 1.В. Качур, В.Ф. Дроботько, Г.А. Городник 

Система монйторингу та морфолопя форми нестащюнарних сигналв внутрипньочерепного тиску 
У статт! розглядаеться розроблена та впроваджена авторами система вимру внутришньочерепного 
тиску (ВЧТ). Пропонуеться для досллдження хвильових форм часових сигналв ВЧТ застосування 
алгоритм1в калмановсько! фальтрацй та калмановського згладжування, що кожного моменту часу 
адаптивно проводять ощнку параметрав авторегресйно! модел1. За допомогою запропонованих алгоритмв 
дослиджено морфологю хвиль ВЧТ за умов гостро! травми головного мозку. 


Р.И. РгоБогКо, 1.Г. Касйиг, Г.Е. ОгоБог’Ко, С.А. СогодшЕ 

Мопкогте Зузет апд Могрво!05у 01 Моп-5айопагу 512па[5 оЁ шегазКиЙ Ргеззиге 

ТБе зу$ет оРшёазКи| ргеззиге (15Р) теазигие е]абогае4 ап4 пар]етещед Бу азот 1$ дезсте4 ш фе 
агасе. Еог фе з4уте оЁТЗР {етрогагу $121а1$ шт йе угауе Ююгил$ И 15 оЁеге4 10 изе а|гогИил$ о Каппап 
ВЯ ганоп ап Каблап зоо те фай ш апу потеп таке адарйуе]у ап аззеззтеп оЁ ашюогестез1юп то4е1 
рагатеегз. 13Р \’ауе тогрпо]озу ипаег зепоиз {гаита оРБгатт 1$ зиеа Бу теап$ оЁ Везе азот. 


Статья поступила в редакцию 09.08.2010. 
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